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Abstract 

The universal quantum work relation connects a functional of an arbitrary observable averaged 
over the forward process to the free energy difference and another functional averaged over the 
time-reversed process. Here, we ask the question if the system is driven out of equilibrium by a 
different Hamiltonian rather than the original one during the forward process and similarly during 
the reversed process then how accurate is the quantum work relation. We present an inequality 
that must be satisfied when the system is driven out by such a trial Hamiltonian. This also answers 
the issue of accuracy of the Jarzynski relation with a trial Hamiltonian. We have shown that the 
correction term can be expressed as the averages of the difference operator between the accurate 
and trial Hamiltonians. This leads to a generalized version of the Bogoliubov inequality for the 
free energy differences. 
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I. INTRODUCTION 



In recent years work relations in nonequilibrium setting have attracted much attention. 
Jarzynski's equality [l| is an important step in that direction, where this provides a relation 
between the distribution of work performed on a classical system subjected to an external 
force and the free energy difference between the initial and final configurations. In particular, 
a finite system is prepared in a state of equilibrium with an environment at temperature T. 
The system is driven with an external parameter from some initial value to a final value. 
During this process the Jarzynski equality connects the equilibrium information about the 
free energy to that of the nonequilibrium processes. The nonequilibrium work relation has 
been proved from a generalized version of the fluctuation theorem [2| . The Jarzynski equality 
has been generalized to quantum domain also HI- Kurchan has considered extension 
of this relation in a quantum setting where he has used measurement based schemes for 
the system of interest at the initial and final time and related the workperformed to the 
energy difference between the final and initial eigenstates [^]. Monnai has derived the 
quantum version of the Jarzynski equality in terms of the microscopic reversibility and has 
explored the equilibrium information about the free energy of a system when it is thermally 
isolated and is driven externally out of equilibrium. The derivation of the work relation by 
Monnai j3] and Tasaki [5] is based on the twice measurements scheme. For more details 
regarding the twice measurements scheme, readers are advised to see @, 0). In addition, 
there is quantum work operator approach to derive the work relation which can give the 
required energy difference when the system is driven out of equilibrium [8|j. The quantum 
analog of the Jarzynski equality has also been derived for closed 0, H, [lOj as well as open 
quantum systems 

Very recently, Andrieux and Gaspard [l2| have proved a universal quantum work relation 
which connects the average of a functional of an arbitrary observable during the forward 
process to that of another functional averaged over the time-reversed process and the free 
energy difference. As shown in that paper many well known results follow from this universal 
work relation. For example, the Jarzynski equality quantum Green-Kubo relations 13]. 
and Casimir-Onsager 14| relations can be obtained from the universal work relation. 

In statistical mechanics approximate methods are very useful when dealing with complex 
systems that consist of many particles. For instance computing the free energy and the par- 
tition function of statistical systems are in general intractable. Only for few special models 
one may be able to calculate such thermodynamic quantities. In these situations one re- 
places the original Hamiltonian (which may be difficult to handle) by a trial Hamiltonian (or 
an approximate Hamiltonian). Under such approximations, proving various thermodynamic 
relations are important and one must try to establish how accurate are those relations. In 
this paper, we ask the question if the system is driven out of equilibrium by a different 
Hamiltonian H'(t) rather than the original one during the forward process and similarly 
during the reversed process then how accurate is the universal quantum work relation. For 
complex systems testing the Jarzynski equality or the Andrieux- Gaspard relation may be 
equally difficult. If we imagine that H'(t) is a trial Hamiltonian then the inequality that 
we are going to prove in the sequel will tell us how accurate is the universal work relation 
or the Jarzynski relation as compared to the actual one. In general, H'(t) may contain dif- 
ferent interaction or several external control parameters. In these cases this inequality will 
be of use in estimating the work performed on the system and the free energy differences. 
In particular, we have shown that the correction term can be explicitly expressed as the 
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averages of the difference operator between the accurate and trial Hamiltonians. Therefore, 
the error in the free energy difference caused by the modelling procedure can be estimated 
by calculating the difference operator. This in turn leads to a generalized version of the 
Bogoliubov inequality for the free energy differences. 

The organization of our paper is as follows. In section II, we introduce the universal 
work relation. In section III, we prove the main inequality for the universal work relation 
under the trial Hamiltonian. In section IV, we address the issue of accuracy of the Jarzynski 
relation under the trial Hamiltonian. Also, we prove a generalized version of the Bogoliubov 
inequality for the free energy differences. In section V, we present yet another useful in- 
equality for the Jarzynski relation in terms of the operator norm. Finally, we conclude with 
some discussions in section VI. 



II. UNIVERSAL WORK RELATION 

Before presenting the main result let us briefly introduce the setting in which the universal 
work relation of Andrieux and Gaspard [l2j was proved. Let us imagine a system whose 
Hamiltonian is given by a Hermitian operator H(t,H), where R can be some external 
parameters that change sign under the time- reversal (such as the magnetic field). Under 
the time reversal operator the Hamiltonian transforms as 0i7(t,R)0 = H(t, — R), with 
being an anti-linear operator and 2 = /. Initially the system is prepared in thermal 
equilibrium with the inverse temperature (3 = where ks being the Boltzmann constant 
and T is the temperature. During the forward process, the system is driven out of equilibrium 
by a time dependent Hamiltonian H(t,R). The initial state of the system is a canonical 
density matrix 



-/3if(0,R) 

p(o) = - m -, (i) 

where Z(0) = Tr[e _ ^^°' R ^] is the partition function, with p(0)^ = p(0),p(0) > and 
Trp(O) = 1. The free energy F(0) of the system at the initial time is given via the partition 
function Z(0) = e~^ F< ^ . We allow the system to evolve from the initial time to some final 
time t = T. The forward evolution is governed by the evolution equation 

ih dUF dt ^ = Hit ' R)f/F(t ' R) - (2) 

If we consider the Heisenberg representation then the observables evolve as A F (t) = 
[/p(t)^A[/i?(t) and similarly for the time-dependent Hamiltonian we have Hp(t) = 
U F (tYH(t,K)U F (t). During the backward process the external parameters are reversed. 
The system is driven by a time-reversed Hamiltonian H(T — t, — R) with a canonical density 
matrix 

Note that the time dependence in the Hamiltonian is rescaled so that the initial time during 
the reverse process corresponds to t — 0. The partition function and the free energy for this 
process are given by the relation Z(T) = Tr[e _/3H( - T ' _:R ^ = e~^ F( ^ T \ The system is allowed 
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to evolve until t = T where the Hamiltonian at the end of the process is H(0, — R)- The 
evolution equation during the backward process is given by 

ih dU *dt R) = HiT ~ *' R)UR(t ' R) - (4) 
Now it can be shown that the forward and the backward evolution operators are related via 

GU F (T - t,R)U F (t,R)Q = U R (t, -R). (5) 

In the Heisenberg picture the observables which have definite parity under time reversal 
obey the following relation during the forward and the backward process 

A F (t) = ±U F (T)QA R (T - t)GU F (T), (6) 

where A R (T — t) = U R (T — t)AU R (T — t). Now, the universal work relation of Andrieux and 
Gaspard states the following. 

Theorem: For an arbitrary time-independent observable A with a definite parity the func- 
tional relation 

( e fi dtX(t)A F (t) e -PH F (T) e mO) )p(olF = e -/3AF (e±/( f d a(T-t)A R (t) )p{T)R (?) 

connects the averages during the forward and the backward processes. Here, X(t) is an arbi- 
trary function, AF = F(T) — F(0) is the free energy difference of the two equilibrium states. 
The average of an operator O during the forward and the backward processes are defined 
via < O > p (o),F— Tr[p(0)O] and < O > p (t),r— Tr[p(T)0], respectively. As mentioned, the 
Jarzynski relation can be obtained from (7). This relation also provides the linear response 
theory of an arbitrary observable. 



III. INEQUALITY FOR QUANTUM WORK RELATION 

In this section, we prove the main inequality for the universal work relation under the 
trial Hamiltonian. Suppose that the system is not driven by the original Hamiltonian but by 
a different Hamiltonian H'(t, R') (we allow the new Hamiltonian to depend on a different set 
of external parameters R'). This can be a trial Hamiltonian or an approximate Hamiltonian 
or it could consist of perturbation to the original Hamiltonian. The question that we are 
concerned here is how stable or how accurate is the universal work relation with such a trial 
Hamiltonian that is different than the original one. To carry out the stability analysis of 
quantum work relation we use adiabatic representation theory where the original Hamilto- 
nian H (t, R) and the trial Hamiltonian H'(t, R') are assumed to have set of instantaneous 
eigenstates at the initial and the final time. Otherwise the results presented here are quite 
general. 

Imagine that the system is driven by a trial Hamiltonian H'(t,TL'). The system starts 
from a canonical density matrix 

-/W(0,R') 

at time t — and evolves as before until time t = T. During the forward process the 
observables evolve according to the trial Hamiltonian as A' F (t) = U' F (tY AU' F (t) and similarly 
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for the time-dependent Hamiltonian we have H' F {t) = U F (t)^H'(t,TL')U' F (t). The backward 
process follows according to the trial Hamiltonian starting from a canonical density matrix 

p -PH'{T,-n') 

AT) = - m -. (9) 

Here, the primed quantities have similar meaning when the system is governed by the trial 
Hamiltonian. The partition function is given by Z'(T) = Tr[e~^ H '^ T ~ R '^]. Now, with this 
trial Hamiltonian one can show that the universal work relation reads as 



{e fT d tX m > F{ t) e -p H > F (T) e PH' { 0) )pi{o)F = e -PAF> {e ±£ dt X { T-t)A> R{ t) )pi{T)R (1Q) 

This connects the functional averages of an arbitrary observable during the forward and 
the backward process when the system is driven by a trial Hamiltonian. Here, X(t) being 
an arbitrary function, AF' = F'(T) — F'(0) being the free energy difference of the two 
equilibrium states. The averages during forward and the backward processes are defined as 
(0) m ,F = Tr[p'(0)O] an d (0) p , (nR = Tr[p'(T)0], respectively. 
What we prove is that the following inequality is satisfied 

^ dt\(t)A' F (t) e -pH F (T) e l3H'(0)^^ F < e -/3AF e /3((y(0)) p(0) , F -(y(T)) p , (T)ii j)^±/ T *A(T-t)A' B (t)) 

'(11) 

where V(t) = H'(t, R') — H(t, R) is the difference between the trial and the actual Hamil- 
tonian. The above inequality can also be stated as 

{e J- ^ m ' A t) e - P H m ^o ))m ^ ^ (e/o T ^ WA - (f) e"^^e^W) p(0) , F c/3((y(0)) ^ w . mT)) ^ p) 

(e ±/ T dtX(T-t)A> R (t)^ /{nR ~ ( e ±£ dtX{T-t)A R (t)} p{T)R 

To prove the above inequality we assume that the time-dependent Hamiltonians H(t) and 
H'{t) have instantaneous eigenstates at time t — and t = T. Let the original Hamiltonian 
satisfies an eigenvalue equation at time t = as if (0, R)|^ n (0)) = E n (0, R)|-0 n (O)). The par- 
tition function at t = can be written as Z(0) = Tr[ e -^ H( - ^} = En(^( )l e ~^ ( °' R) l^n(0)>- 
During the forward process consider the diagonal elements of the canonical density matrices 
p(0) and p'(0) in the eigenbasis \if) n (0)). Thus, we have two probability distributions p n (0) 
and p' n (0) as given by 



Pn(0) 



e -/3E„(0,R.) 

W) 



-f}H'(0,R') 

P'n(0) = (^n(O)l - Z/(Q) \M0))- (12) 

For these two probability distributions p„(0) and p' n {0) the relative entropy of p with 
respect to p' (also called the Kullback-Leibler distance) is defined by 

R(p,p') = J2Pn(0)hg^- ) >0. (13) 
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It is a convex function of p n (0) and is always non- negative and equals zero only if p n (0) = 
p' n (0). Using this, we can write the inequality as 

]Tp n (0) log Pn (0) > ^p„(0) logK(O). (14) 

n n 

Now, using the expressions for p n (0) and p' n (0), we can write the above inequality as 
£p n (0)[-/3£ n (0,R)-logZ(0)] >J2Pn(0)(Mm-PH'(0,R')-logZ'(0))\i; n (0)).(15) 

n n 

In the above equation, on the right hand side we have used the Jensen inequality, i.e., 
(e°) > which gives us 

logK(O) = log^nWK-^'^RO-logZ^O))!^^))] 

> (^(0)|(-/?/J'(0,R0-logZ'(0))|^ n (0)) (16) 

Therefore, we have the following inequality that gives a lower bound on the partition function 
with the trial Hamiltonian at the initial time of forward process 

Z'(0) > Z(0)e _/3Tr[p(0) (tf'(°.R')--ff(o,R)]_ (17) 

Next we derive an upper bound on the partition function with the trial Hamiltonian 
at the initial time of the backward process. During the backward process with a trial 
Hamiltonian the system starts from the density matrix as given in (9). Let us assume 
that the trial Hamiltonian satisfies an instantaneous eigenvalue equation at time t = T, i.e., 
H'(T, -R')K(T)) = K( T , -R')K( T ))> where K( T , -R-') is the eigenvalue. Now consider 
the diagonal elements of the canonical density matrices corresponding to the Hamiltonian 
H(T, — R) and H'(T, — R') in the eigenbasis of \tp' n (T)). These probability distributions are 
given by 

-I3H(T-R,) 

Pn(T) = mn z(T) w n {T)) 

-/3if'(T-R') f> —[3E' n (T,—~R!) 

Pn(T) = «(T)| zf{T) K(T)) = ZI[T) . (18) 
Let us consider the relative entropy of p' n (T) with respect to p n {T) which is given by 

? K(T)1 %ifi aa (19) 

Again using the non-negative property of the relative entropy we can write the above in- 
equality as 

J2Pn(T) logpUT) > Y.P'n( T ) logPn(T). (20) 

n n 

Using the expressions for p' n (T), p n (T) and the Jensen inequality as before, we have 

Y,Pn(T)[-PK(T, -R 1 )- log Z'(T)]> 

n 

E Pn( T ) MT) I (—/3H(T, -R) - log Z(T)M(T)> . (21) 
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This leads to an upper bound for the partition function with the trial Hamiltonian during 
the time-reversed process as given by 



Z\T) < Z(T)e-W™ h'(t,~r')-h(t,-K) )]_ (22) 

From these two inequalities (17) and (22) we have 

Z '( r ) < Z ( T h mv(Q)) n cnvp-<v(T)) r , l ^ ^ (o*\ 

Z'(0) ~ Z(0) ' 1 ' 

where (V{0)) p{0)>F = Tr[p(0)( H'(0,R') H(0,R) )] and (V(T)) p , (nR = 

Tr[p'(T)( H'(T,-R') - H(T, — R) )]. By noting the fact that = e~^ AF ' and 

: |^y = e~ l3AF , we have the main inequality for the universal quantum work relation with a 
trial Hamiltonian as given in Eq.(ll). Hence, the proof. 



IV. JARZYNSKI RELATION WITH TRIAL HAMILTONIAN 

In this section, we address the question of the accuracy of the Jarzynski relation when the 
system is driven away from equilibrium by a trial Hamiltonian. Interestingly, the quantum 
version of the Jarzynski equality follows from the universal work relation of Andrieux and 
Gaspard [l2|. If we set A = in Eq.(7), we obtain 

{e -m F {T) e mo) )F = {e -m ) = c -/3af (24) 

where the left hand side represents the average of the exponential of work performed on 
the system during the forward process and AF is the equilibrium free energy difference 
F(T) — F(0). Now if one asks the question what happens to the Jarzynksi equality if we 
drive the system with a different Hamiltonian H'(t) rather than H(t). Here, H'(t) could be 
a trial Hamiltonian, or could be some perturbation. Under such a situation how stable is 
the Jarzynski relation or the nonequilibrium work relation. Our inequality can answer this 
question. It is clear that if we drive the system with a different Hamiltonian H'(t, R') then 
we will have the Jarzynski equality 

( e -m'AT) e W'{*)) F = ( e -pw') = e -/^. (25) 

From our inequality (11), if we set A = 0, we will obtain the following inequality for the 
Jarzynski relation 

, e -/3B' F (T) e l3H'(0)\ = e -pAF' < e -/3AF e /3[{y(0)> p(0) , F -{V(T)) p , {T) , fl ] 

or {e- m> ) < ( e -P w } e P [{vm »^ F - {v{T)) ?'m>K ] . (26) 

This tells us that the work performed by a trial Hamiltonian in going from the initial 
configuration to the final configuration must respect this bound. Eq.(26) is important in its 
own right. This shows that the correction term can be explicitly expressed as the averages 
of the difference operator between the accurate and trial Hamiltonians. By calculating the 
quantity [(V"(0)) p (o),f — (^ / (7 1 ))p'(r),i?] the error in the free energy difference caused by the 
modelling procedure can be estimated directly. 
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Furthermore, using the inequality (26), we can prove a generalized version of the Bogoli- 
ubov inequality. One version of the usual Bogoliubov inequality [l5[ can be stated as follows: 
Let H be the original Hamiltonian and let H' be the trial Hamiltonian with the condition 
that tr(p'H') = tr(p'H), i.e., the averages of both the Hamiltonians in the canonical state p' 
are same, then F < F'. Here, F is the free energy of the original Hamiltonian and F' is the 
free energy of the trial Hamiltonian. The Bogoliubov inequality has important application 
in the mean field theory. For example, by a variational method if we can minimize the free 
energy of the trial Hamiltonian we can get a better approximation to the exact free energy. 

To prove the generalized version of the Bogoliubov inequality, we note from (26) that 

AF' - AF > [(V(T)) p , (nR - (V(0)) p(0) , F }. (27) 

Now, if we impose the condition that H(t) and H'(t) have same averages in the original and 
trial canonical state at time t = and t = T, then (V(T)) p >(t),r = and (V"(0)) p (o),f — 0. 
Then it follows that 

AF < AF'. (28) 

This can be regarded as time-dependent generalization of the Bogoliubov inequality which 
bounds the change in the free energy of original Hamiltonian with respect to perturbed one. 
This is an another important result. This inequality not only holds close to equilibrium but 
also hold in those situations which are far-from-equilibrium regime. 

V. INEQUALITY WITH OPERATOR NORM 

In this section, we prove an inequality for the Jarzynski relation which involves the norm 
of the operator V(t). This may be useful when we do not have to calculate the averages of 
the operator V(0) and V(T) during the forward and backward processes, corresponding to 
the actual and trial Hamiltonians. 

Consider the partition function for the trial Hamiltonian during the initial time of the 
reverse process. We have 

Z'(T) = Tr[e-^'( T >- R,) ] < Tr[e^ T ^e^ v ^}, (29) 

where we have used the inequality Trfe^" 8 ] < Tr[e A e B ] for all self-adjoint operators. 
On expressing the trace using the eigenbasis of H(T, — R) we have 

Z'(T) < ^(i(T)| e -W- R )e-^|i(T)) 

n 

= ^(^(F)|e-^- R )|^(F))(^(T)|e-^ T )|^(T)). (30) 

n 

Now, we use the inequality \\e A \\ < e"" 4 " which holds for all operators. The operator norm 
of a linear operator A is defined as \\A\\ = sup(||^|i_ 1 \||-A(^>)||. (Note that in the matrix 
notation, A(ip) = A\ip) and \ \A\ \ is equal to the square root of the largest eigenvalue of the 
symmetric matrix A T A .) 
Using this we have 

Z'{T) < ^(^(T)|e-^^- R )|^(T))e^l^)H 

n 

= Z(T)ef>" v <™. (31) 
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This gives an upper bound for the partition function for the trial Hamiltonian at the initial 
time of the reverse process. 

Now, let us consider the partition function Z(0) with the Hamiltonian H(0) during the 
initial time of the forward process. We have 

Z(0) = Tr[e-^°< R )] < Tr[e^ H '^e^}, (32) 

where we have again used the inequality Tr[e A+jB ] < Trfe^e 5 ] which holds for all self-adjoint 
operators. We express the trace using the eigenbasis of H'(0,H). This gives us 

z(o) < E«(°)l e "^' (0 ' R ' )e/3y(0) l<(°)) 

n 

= E«( )l e "^' (O ' R ' ) l<(°))^(0)|e^^|^(0)). (33) 

n 

Now, we use the inequality ||e A || < e^ A ^ which holds for all operators. Using this we have 

Z(0) < E(^(0)|e-^'(°' R '^;(0))e^l^)ll 

n 

= Z'(0)e /3 H y (°)". (34) 

This gives a lower bound for the partition function with the trial Hamiltonian during the 
initial time of the forward process. 

On multiplying two inequalities (31) and (34), we have the new inequality for the Jarzyn- 
ski relation under the trial Hamiltonian as given by 

(e~ pw ') = e _/3AF ' < e ~ fiAF e mm e m( - T)l1 
or (e-W) < ( e -^) e /3||v(o)|| e /3||y ( T)||_ (35) 

This inequality may be often useful where we do not have to calculate explicitly the averages 
of the operator V for the actual and trial states. Eq. (35) provides the exact estimation of 
the free energy difference between the initial and final states when the actual Hamiltonian 
is not the same as the trial or approximate Hamiltonian. With the help of the free energy, 
all the quantum equilibrium properties like different phases of a system and their physical 
properties like the phase transitions etc. can be exactly calculated. This relation is expressed 
in terms of an inequality which involves a correction term when the actual Hamiltonian is 
not same as the trial or approximate Hamiltonian. The correction term is expressed in the 
form of the norm of the operator V at the initial and final time. In some physical situations 
it may be difficult to calculate exactly the averages of V(T) or V(0). Also, the free energy 
of a quantum system relative to an arbitrary reference state is often difficult to determine. 
In such situations, Eq.(35) may provide a practical tool to extract the quantum equilibrium 
information even if the averages of V is explicitly not known. The norm of V at the initial 
and final time is enough to give the required estimate. 

VI. CONCLUSIONS 

Testing various thermodynamic relations for complex systems are not always amenable. 
One has to model the actual Hamiltonian of the system by a trial Hamiltonian and establish 
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how accurate are those relations. In this paper, we have investigated the stability of the 
universal quantum work relation of Andrieux and Gaspard under a trial Hamiltonian and 
proved an inequality that must be respected if the system is driven away from equilibrium 
by a trial Hamiltonian rather than the actual one. In turn, our inequality also tells us 
how accurate is the famous Jarzynski relation under such trial Hamiltonians. We have 
shown, in particular, that the correction term can be explicitly expressed by the averages of 
the difference operator V between the accurate and the trial Hamiltonians. Therefore, the 
free energy difference caused by the modelling procedure can be estimated by calculating 
the difference operator. We have also obtained a generalized version of the Bogoliubov 
inequality for the time- dependent trial Hamiltonians. This tells us that the change in the 
free energy corresponding to the actual Hamiltonian is always less than or equal to the 
change in the free energy in the trial case. Further, we have given an inequality using the 
operator norm where the correction term can be expressed in the form of the norm of the 
operator V at the initial and final time. In future, we plan to apply this inequality to some 
realistic experiments. We hope that this inequality will be very useful in testing quantum 
work relations for many body systems and will open up further thoughts in the area of 
quantum work relations and fluctuation theorems. Also, this inequality may be applied to 
probe quantum nonequilibrium phenomena in glassy systems and many body systems. 
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